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We calculate the distribution of eigenfunction amplitudes 
and the variance of the "inverse participation ratio" (IPR) in 
disordered metallic samples. A relation is established between 
the distribution function of IPR and that of "level velocities" 
- parametric derivatives of energy levels with respect to a 
random perturbation. 

PACS numbers: 71. 55. J, 05.45. +b 



I. INTRODUCTION. 

Statistical properties of disordered metals have at- 
tracted a considerable research interest last years. It was 
understood that the old problem of a quantum parti- 
cle moving in a quenched random potential considered 
earlier in the context of Anderson localization and meso- 
scopic phenomena Q exemplified a particular class of 
chaotic quantum systems and had much in common with 
such paradigmatic problems in the domain of Quantum 
Chaos as quantum billiards 0. The Wigner-Dyson en- 
ergy level statistics first found in the framework of ran- 
dom matrix theory (RMT) || and considered to be a 
"fingerprint" of quantum chaotic systems 0] was shown 
to be relevant for disordered metals as well This 
fact gave a boost to a broad application of RMT results 
for qualitative and quantitative description of mesoscopic 
conductors and stimulated a common interest to statis- 
tical characteristics of spectra of disordered systems [|| . 

At the same time less attention was paid to statisti- 
cal properties of eigenfunctions in disordered or chaotic 
quantum systems. Recently, however, the distribution of 
eigenfunction amplitudes was shown to be relevant for de- 
scription of fluctuations of tunnelling conductance across 
the "quantum dots" j^] as well as for understanding some 
properties of atomic spectra (To). Besides, a so called 
"microwave cavity" technique has emerged |HJ as a lab- 
oratory tool to simulate a disordered quantum system. 
This technique allows to observe directly eigenfunctions 
spatial fluctuations and was used in to study experi- 
mentally the eigenfunction statistics in weak localization 
regime. All these facts make the issue of eigenfunction 
statistics to be of special interest and are calling for a 
detailed theoretical consideration. 

In order to characterize eigenfunction statistics quan- 
titatively, it is convenient to introduce a set of moments 



I a = J \ijj(r)\ 2q d d r of eigenfunction local intensity |V , ( r )| 2 
|13| . The second moment 1% is known as the inverse par- 
ticipation ratio (IPR). This quantity is a useful measure 
of eigenfunction localization: it is inversly proportional 
to a volume of a part of a system which contributes ef- 
fectively to eigenstate normalization. For completely "er- 
godic" eigenfunctions covering randomly, but uniformly 
the whole sample I2 oc 1/V, with V being a system vol- 
ume. If, in contrast, the eigenfunctions are localized, 
i.e. concentrated in a region of linear size £, the mean 
IPR scales as I2 oc £~ D2 where D 2 is an effective dimen- 
sion which can be different from a spatial dimensionality 
d because of a multifractal structure of eigenfunctions 
p3[ . Correspondingly, IPR fluctuations reflect level-to- 
level variations of eigenfunction spatial structure. 

The most complete analytical study of statistical char- 
acteristics of eigenfunctions was performed for the cases 
of Od system s |14| , |26|| , as well as for strictly Id 15 and 



quasi Id JT^,|l8| geometry Some analytical results were 
obtained also for a system in the vicinity of localization 
transition in the dimensionality d = 2 + e, e 1 |i~5| ] 
as well as for d — > 00 (T^|. Let us note that in Ref. 
Jl^ , p7| - [l9|j2^ | the supersymmetry method was used which 
is a very powerful tool to study distribution functions 
of various quantities characterizing eigenfunctions statis- 
tics. 

In the present article we address systematically the is- 
sue of the eigenfunction statistics for arbitrary spatial 
dimensionality d in the weak localization domain. In 
the leading approximation (which neglects spatial struc- 
ture of the system and treats it as a zero-dimensional 
one) these statistics are described by the RMT which 
predicts a Gaussian distribution of eigenfunctions am- 
plitudes ip( r ) ]l4l|lS| ] for systems with unbroken or com- 
pletely broken time reversal symmetry |pof . It is known 
since the paper by Altshuler and Shklovskii JtJ that the 
diffusion motion of a particle in a metallic sample pro- 
duces deviations of spectral statistics from that predicted 
by RMT. To our best knowledge, the analogous problem 
for the eigenfunctions statistics for 2D or 3D systems has 
never been studied. It is just considered in the present 
paper. We use a recently developed method |2l[] which 
is based on the supersymmetry technics |^,|(| and com- 
bines a perturbative elimination of fast diffusive modes 
(in spirit of renormalization group ideas) and consequent 
non-perturbative evaluation of a resulting Od integral. In 
this way, we calculate the deviations from the Gaussian 
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distribution of ip{r) in mesoscopic metallic samples. We 
calculate also the variance of the IPR which turns out to 
be of order of 1/g 2 , where g is the dimensionless (mea- 
sured in units of e 2 /h) conductance of the sample. In 
the final part of the paper we establish a relation be- 
tween the IPR distribution Viil-i) and the distribution 
V v (v) of the so-called "level velocities" v, which charac- 
terize a motion of energy levels as a response to a small 
random perturbation. It turns out that the finite value of 
IPR variance gives rise to deviations of the form of V v (v) 
from the Gaussian. The latter form is known to be typ- 
ical for completely "ergodic" systems with unbroken or 
completely broken time reversal invariance p2|-p4|] . 



II. DISTRIBUTION OF EIGENFUNCTION 
AMPLITUDES AND IPR FLUCTUATIONS. 

In order to calculate the distribution of eigenfunction 
amplitude and to find the IPR variance we use the fact 
that relevant quantities can be expressed in terms of cor- 
relation functions of certain supermatrix a- model §,§. 
A q uite general exposition of the method can be found in 
p5| and is not repeated here. Depending on whether the 
time reversal and spin rotation symmetries are broken 
or not, one of three different cr-models is relevant, with 
orthogonal, unitary or symplectic symmetry group. We 
consider first the case of the unitary symmetry (broken 
time reversal invariance) and then present the results for 
two other cases. 



A. Unitary symmetry. 

The expressions for I q and I 2 (the bar standing for dis- 
order averaging) in terms of the cr-model read as follows 
0HI1|: 
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I? = — lim e 3 -—r 

2 6V e^O du 2 



f DQexp{-fM(u,Q)} ; (1) 

u=0 J 

[ DQcxp{-T (2 \u,Q)} ; (2) 



where 



<r'rt -Str(VQ)' 



+eStr(AQ) + uStr«(QAfc)} 
Q = T _1 AT ; A = diag(l, 1,-1,-1) ; 



k = diag(l, — 1, 1, —1) ; 



- = ttvV/A 



(3) 



Here T is 4 x 4 supermatrix belonging to the coset space 
U(l,l | 2), T> is the classical diffusion constant, v is the 
density of states and V is the system volume. 

Generally speaking, the RMT predictions are applica- 
ble to a disordered metallic system under the following 



conditions: L ^> I; E c ^> A, where L is the system size, 
I is mean free path , E c = KD/L 2 is the Thouless en- 
ergy and A is the mean level spacing. Eigenfunctions for 
such systems are known to be ergodic with amplitudes 
ip(r) being uncorrelated ( for | r — r' |> I ) Gaussian dis- 
tributed complex (real) variables for broken (unbroken) 
time-reversal symmetry respectively. This immediately 



gives 4 = ql/V^ 1 and I 2 = I q , where the 



superscript u refers to the unitary symmetry. 

In the framework of the a— model formalism these re- 
sults can be easily reproduced if one neglects any spatial 
variation of the supermatrix field Q(r). Then eqs.(]l|), (Q) 
are reduced to integrals over a single supermatrix which 
can be evaluted exactly. The corrections to RMT results 
have a form of a regular expansion in small parameter 



A/E c 



ff 



The systematic way to construct such an 



expansion can be briefly outlined as follows pij| . The 
matrix Q(r) is decomposed as 

Q{r) = Tz l Q{r)T Q , (4) 

where To is a spatially uniform matrix and Q describes 
all modes with non-zero momenta. When A -C E c , the 
matrix Q fluctuates only weakly around the value Q = A. 
It can be parametrized as H 



Q = A(1 
= A ( I 



W/2)(l - W/2)- 1 
W 2 W 3 



w 



w 4 



(5) 



where W is a block off-diagonal supermatrix representing 
independent fluctuating degrees of freedom. Substituting 
this expansion into eqs.([l])-(|3|) and integrating out the 
"fast" modes one obtains an expression for renormalized 
functional !F^^(u, Qo), where Qq = T _1 ATo is an i — 
independent matrix (zero mode): 

e -^ f (u,Q ) = i e -^"Hu,Q ,W)+lnJ(W)\ w _ (g) 

Here J(W) is the Jacobian of the transformation (||), (||) 
from the variable Q to {Qo, W}; and (. . .)\y denotes 
integration over W. Now one expands the exponential in 
(||) and calculates integrals over the fast modes W using 
the Wick theorem and the contraction rules : 

(StrW(r)PW(r')R) 

= n(r, r'XStrPStri? - StrPAStriJA) ; 
(Str[W(r)P]Str[W(r')R]) 

= IL(r,r')StT(PR- PARA); (7) 

where P and R are arbitrary supermatrices. Here 
H(ri, T2) is the diffuson propagator; for an isolated sam- 
ple this propagator has the following form: 



n(r l7 r 2 ) 



U(q) = 



E 



cos(qri) cos(qr 2 ) U(q) 
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2ttvV Vq 2 + e 
m = 0,±l,±2,... ; 



m n d 

q = 7T — 



L 



(8) 
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where the system is thought to be of the size L\ x L 2 x 
... x Ld- Finally, the integrals over Qo are performed 
exactly, using the technique developed in ||[| . 
Applying this method to eqs.(0), (^) one obtains: 
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(u) _ gi 



1/9-1 
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(9) 



where <? = 2-KvT)L d 2 is the conductance of the sample. 
The value of the coefficient ai, 



ai 



<E n (< 



i 



E 



7r^ — ' nr 

ri;>0: n 2 >0 1 



(10) 



depends on the spatial dimension and is equal to a\ = 1 /6 
in quasi lei systems. For d > 2 the corresponding sum 
over momenta q diverges at large \q\ and is to be cut off 
at |q| ~ l^ 1 . This gives a\ = ^-\sxL/l for d = 2 and 
a\ <x L/l for d = 3. 

Knowing all the moments iq it is easy task to restore 
the whole probability distribution V{y) of the eigenfunc- 
tion local intensity y = V\ip(r)\ 2 : 



T {u) (y) = e- y 



l + -(2-4y + y 2 ) 
9 
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The leading terms here reproduce the well-known RMT 
result JU; the rest is the weak localization correction. 

Now we turn to the consideration of IPR fluctuations. 
It turns out that IPR variance is of order of l/g 2 . Thus 
the expression (j^) is insufficient for our needs and should 
be extended to the next order. By the same method we 
hnd: 
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(12) 



(13) 



Here the coefficient a 2 is dchncd as a 2 = 5 2 X) q n 2 (q) 
and is equal to 



1 

«2 = —7 



1 



n 4 >0; n 2 >0 



2^,2 



(14) 



the sum being convergent for d < 4. In particular, for 
d = 1,2,3 we have a ± = 1/90 ~ 0.0111, a 2 ~ 0.0266 and 
03 ~ 0.0527 respectively. 

Thus, we find the following expression for the relative 
variance of the IPR distribution: 



5 (u) (7 2 ) 



rO) 



r(«) 



r(«) 



5a 2 
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(15) 



B. Orthogonal and symplectic symmetry. 

For the systems with unbroken time reversal symme- 
try (orthogonal and symplectis cr-models) the calculation 
goes along the same lines. The main difference is in the 
form of the contraction rules which in this case include 
an additional (Cooperon) contribution ||^9|,[50| and can 
be presented in the following form: 

(StrW{r)PW(r')R) 

= 2n(r,r')[StrPStri? - StrPAStriiA 

+aStr(PARA - PR)} ; 
(Str[W(r)P]Str[W(r')R]) 

= 2H(r, r')Str(PR - PARA -PR + PARA) (16) 

where a = +1 (—1) for the orthogonal (resp. symplec- 
tic) symmetry, and R = M T R T M , with M being a su- 
permatrix representing the time reversal transformation 
and satisfying M T M — 1, M 2 = ak (see [[30| for more 
details). 

It turns out however that the formulae (||), (|l5| ) can 
be generalized to the orthogonal and symplectic cases 
without additional calculations. The leading weak- 
localization correction to I q originates from the one- 
diffuson-loop diagram which is proportional to q(q — 
l)a\/g. The numerical coefficient can be restored by us- 
ing the results we obtained earlier for the case of quasi- 
1D systems 025): 



r(lB) 



= I^{l + lq(q-l)x 

+ j^g(g - i)(g - 2)(3« - i)x 2 + o(x 3 )} , (17) 

where x — 2/((3g); (3 = 1(2,4) for the orthogonal (resp. 

unitary, symplectic) symmetry, and Iq GE%> denotes the 
value of I q for the corresponding Gaussian ensemble. We 
see that the quasi- ID result is in agreement with our 
general formulae (||), ([l2]). We restore therefore the coef- 
ficients of the leading corrections to I q in the orthogonal 
and symplectic cases: 



«{ 1+ | 9(9 _ 1)+ o(i)} 



Asp) _ (g + 1)! J, a 



(18) 
(19) 



The corresponding corrections to distributions of eigen- 
functions intensities are as follows lETT]: 



,-y/2 



V (o \y) — 
s/2-ny 

V [sv \y) = Aye~ 2y 



1 + — (3 - 6y + 2y 2 ) 
9 



(20) 
(21) 



As to the IPR variance, 5(i r 2 ), the leading contribution 
to it is given by the diagrams with a two-diffusons-loop 
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and is therefore proportional to a 2 / g 2 . The coefficient 
can be again restored by comparison with the quasi-lD 
result dUl]: 



5 {1D) {h) 
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45/3 : 



2„2 



(22) 



In the quasi-lD case a 2 = 1/90, so that ( |22| ) agrees with 
( |l5| ) . For orthogonal and symplectic systems of a general 
dimensionality, we find 



<5 (sp) (/ 2 ) 



32a 2 / 5 2 
2a 2 / ff 2 



(23) 
(24) 



We have checked the above results by direct calculation 
in the orthogonal case. It should be noted that in order 
to get the correct results for 7 2 and I 2 to l/<? 2 order 
in this case, one should take into account the Jacobian 
J(W) in eq.(||) (whereas in the unitary case this Jacobian 
contributes only starting from the l/g 3 order). For the 
orthogonal symmetry, we find 



LnJ(W) 



StvW 2 + 0(W 4 ) 



Carrying out the calculations we obtain 



r(°) 



3 
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l^ + l(12 a 2 -30a 2 ) + O(l 
9 9 \9 6 
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-\{40a 2 -28a 2 ) + O (\ 



(25) 



(26) 



(27) 



and consequently, S^'fa) = 32a 2 /<? 2 , in agreement with 



C. Discussion. 

Let us now give some comments on the results obtained 
in this section. Eqs.([ll]), (pp|), ( pl| ) present leading (one- 
loop) corrections to the distribution of eigenfunction am- 
plitudes which are due to the weak-localization effects. 
It is easy to see that the corrections lead to an increase in 
probability to find a value of the amplitude considerably 
smaller or considerably larger than the average value. It 
is exactly what one would intuitively expect assuming a 
tendency of cigenfunctions to localization. It is interest- 
ing to note that the correction is of the same type for all 
three ensemble, in contrast to the well known difference 
in weak-localization corrections to the conductivity. 

Deviation of the eigenfunctions amplitude distribution 
from its RMT form was observed very recently in an 
experiment with a disordered 2D microwave cavity [ jl2| . 
The form of the correction reported in is in agree- 
ment with our formula (p0|). 



It should be noted that Eqs.Q, @, @ are valid 
for not too large values of y: y < y/g/ax, when the cor- 
rection term is small as compared to the leading one. For 
larger y (i.e. in the far "tail"), the distribution V{\) dif- 
fers strongly from its RMT form and can not be found 
by the method used in this paper. At present we are 
able to present the result for the "tail" of V(y) for the 
quasi-lD geometry only. This can be done using the 
exact analytical expressions for V(y) we have found ear- 
lier for quasi-ld systems (see Eqs.(16), (17) of Ref. JlJ] 
or Eqs.(79), (89) of Ref. ||). We find three different 
regions of the variable y: 



P<%)~e-» l + -(2-4y + ? y 2 ) + 
L 6 



; y < x- 1 ' 2 (28) 



V {u \y) ~ exp 



-y H y x 



V2 < < 



y^x' 1 (29) 



~ exp - / l\fy/x ; y > x^ 1 



(30) 



where x = l/g. For completeness, we repeated again the 
result of preceding subsections, Eq. (p8|) , which presents 
a relatively small correction to the RMT formula. In 
the intermediate region j29|), the correction in the expo- 
nent is small compared to the leading term but much 
larger than unity, so that V(y) ^> Vrmt{v) though 
lnV(y) ~ hi'pRMT(y)- Finally, in the region (|30|), the 
distribution function V(y) has nothing to do with its 
RMT form. To understand the physical origin of the 
asymptotical behavior (|3(]) we note that it has exactly 
the same form as in the case of a system of an infi- 
nite length, when all states are localized (see Eq.(21) 
of Ref. @ or Eq.(114) of Ref. §|). In the quasi-lD 
case x = l/g = L/L c , where L is the sample length 
and L c is the localization length in an infinite system. 
Thus, the ratio y/x in the exponent in Eq. (p0|) is equal 
to y/x = V c \rp(r)\ , where V c is the localization volume. 
Therefore, the asymptotics ( |30| ) is controlled by the cor- 
relation length L c and describes a probability to have a 
quasi-localized "bump" with an extent much less than L c 
p8|. We stress again that the formulae (p^), ( j30| ) have 
been derived for quasi-lD systems; the problem of their 
generalization to higher d is still open. 

The obtained value of the IPR fluctuations £(/ 2 ) ~ 
l/g 2 oc L 4 ~ 2d is much larger than the RMT result 
<5(/ 2 ) cx 1/V ~ L~ d . These fluctuations have much 
in common with the universal conductance fluctuations 
which have the same relative magnitude. It is interest- 
ing to note that extrapolating this result to the Ander- 
son transition point, where g ~ 1, we get £(/ 2 ) ~ 1, 
so that the magnitude of fluctuations is of order of the 
mean value. Thus, despite the fact that the wave func- 
tion is extended through the whole sample, there is no 
self-averaging of the IPR at the mobility edge. 
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III. DISTRIBUTION OF LEVEL VELOCITIES. 

Let us now demonstrate that the IPR fluctuations are 
intimately related to level response characteristics of dis- 
ordered system subject to a random perturbation. 

The issue of the energy level motion under an ex- 
ternal perturbation was under intensive study recently 
p2| , p3 31 3^]. In the paper j2^| the statistical properties 
of the "level velocities" (LV) defined as v n — dE n /da, E n 
standing for a given energy level and a being some con- 
trol parameter, were studied for systems with completely 
"ergodic" eigenstates. In particular, the LV distribution 
V(v) was shown to be Gaussian within the zero- mode ap- 
proximation for supersymmetric er— model. On the other 
hand, in the recent paper J32| the LV distribution was 
calculated analytically for a long quasi Id system sub- 
ject to a random perturbation. It was found that in 
the limit g — > which corresponds to complete localiza- 
tion of eigenstates, V(v) is essentially nongaussian. The 
mean squared LV (v^) was found to be of order of in 



The explicit form of fed(r) for d = 2, 3 can be found in 
p3| ; for our purposes it is enough to mention its asymp- 
totic behaviour: 



fc rf (r) ~ 1 
k d (r) < 1 



fc/r < 1 
fc/r > 1 



(34) 



where kf is the Fermi momentum. 



Substituting ( p3| ) in (32) and using (|34| ) we find 

(^2) = Wa I 2 ; Wo = c f d d rW (r) , (35) 



where I 2 is the mean IPR introduced above and 
c = 1 (1/2) if k f ( < 1 (resp. fe/C > 1^ 

Therefore, the mean squared LV (v%) is proportional 
to the mean IPR I 2 . In the infinite quasi Id system I 2 
is inversly proportional to the localization length £ in 
agreement with the behaviour (v^) oc found in |32| . 

This consideration can be extended to higher LV mo- 



verse localization length £-\ whereas the RMT predicts ments ( w » > as welL After some manipulations we find 



In order to get understanding of these facts in broader 
context let us consider a disordered system as described 
by a one-particle random Hamiltonian TL. Let us impose 
on it some random perturbation W so that perturbed 
Hamiltonian is TL + aW, with the parameter a control- 
ling the strength of perturbation. We assume that both 
the Hamiltonian TL and the perturbation W belong to 
the same pure symmetry class [p4|. In the sequel, we 
denote the averaging over the realizations of W by (...) , 
which should be distinguished from the averaging over 
unperturbed Hamiltonian TL denoted by bar. The level 
velocity v n corresponding to an energy level E n can be 
found within the conventional perturbation theory as 



the general relation between the two distributions: 



V v (v) = (S(v- 



v n )) 
dl 2 



[2ttw I 2 } 1/2 



exp 



2w I 2 



Vi{h) (36) 



d d r / d d r'Wrr>r n (r)Mr') 



(31) 



where W rr ' are matrix elements of the perturbation. 
To be specific, these matrix elements are supposed to 
be independent Gaussian distributed random variables 
with the mean value equal to zero and the variance 
(W* r ,W r r>} = W {\r - r'\). We also make a natural 
assumption that the perturbation W is of finite range 
i.e. Wo (1^1 ^ C) is assumed to vanish. 
By using eq.(pl|) one easily obtains: 



The following comment is appropriate here. For a 
given realization of the random perturbation W one can 
define the mean level velocity tvT (note that the averag- 
ing here goes only over the unperturbed Hamiltonuan TL, 
but not over W). It follows from eq.(|3^) that de- 
fined in such a way shows Gaussian fluctuations around 
zero from one realization of the perturbation to an- 
other. The mean square (u^ 2 ) can be easily found to be 
(v 2 ) = V^ 1 J d d rW (r)\G{r)\ 2 and is proportional to the 
inverse system volume V -1 . At the same time the distri- 
bution function V u (u) of the quantity u — v — turns 
out to be self-averaging, i.e. its form is the same for any 
realization of the perturbation W. The function V u (u) 
can be shown to satisfy the relation similar to eq.(pq): 



V u {u) 



dh Vi{h) 



[2n(w I 2 - (v 2 ))] 



exp 



~2(w h- (v 2 })_ 
(37) 



d d r / d d r'W (\r 



Mr)Mr') I 2 (32) 



To simplify eq.(|3^) we need to know the short-scale be- 
haviour of the correlator |'0^(r)-0 2 l (r / )|. We have EjjJ: 



\^{rW{r')\ = - [1 + k d (\r r'|)] |^(r)| , (33) 



where kd (r) = ImG(r) /(tt^) 2 and G(\r — r'\) = (r\(E ■ 



ny 



is the Green function of the Hamiltonian TL. 



The relation (|36|) is then recovered as a convolution of 
eq.(|37]) with the Gaussian distribution of v n . 

Eq.(^6|) shows that the Gaussian form of V v (v) found 
earlier for generic chaotic systems |^,^3| is a consequence 
of the IPR being a selfaveraging quantity within zero- 
mode approximation to nonlinear a— model p4j . For 
disordered systems with finite conductance g the func- 
tion V v (v) becomes nongaussian due to mesoscopic fluc- 
tuations of cigenfunctions. As a natural measure of the 
deviations from the Gaussian form one can use the cu- 
mulant: 
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(IP =3^ 5 - = 3 TO « ? (38) 

With increasing disorder the conductance decreases 
and the system passes to the strong localization regime 
g -C 1. For this regime the explicit expression for both 
distribution functions Viil-i) and V v (v) were found in 
quasi Id systems in the course of independent calcula- 
tions (l7H|] (see als o PI ) and ||. The results indeed 
satisfy the relation (p6|). 

It could be interesting to study the functions V v (v) and 
Viil-i) also in d = 2 and in the vicinity of the Anderson 
transition for d > 2 where their form is expected to be 
universal. Relation ( ^6|) may provide then a convenient 
background for investigations of the statistical properties 
of eigenf unctions. In particular, the behaviour of the 
mean square level velocity (v 2 ) in the critical region can 
be used for extracting a value of the anomalous effective 
dimension Z?2- 

IV. SUMMARY 

In this paper we have studied deviations of the eigen- 
function statistical characteristics in a disordered metal- 
lic sample from those predicted within the Random Ma- 
trix Theory. The found correction to the distribution of 
eigenfunction amplitudes reflect the influence of weak lo- 
calization to an eigenfunction. The relative magnitude 
of fluctuations of the inverse participation ratio is of the 
order of 1/g 2 , where g is the dimensionless conductance, 
in close analogy with the conductance fluctuations. We 
have also revealed the relationship between the meso- 
scopic fluctuations of the inverse participation ratio and 
the form of the "level velocity" distribution in a system. 
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